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ABSTRACT 

The role of the acoustic continuum associated with compact sources in the Sun's interior wave field 
is explored for a simple polytropic model. The continuum produces a near-field acoustic structure 
- the so-called 'acoustic jacket' - that cannot be represented by a superposition of discrete normal 
modes. Particular attention is paid to monochromatic point sources of various frequency and depth, 
and to the surface velocity power that results, both in the discrete /- and p-mode spectrum and in the 
continuum. It is shown that a major effect of the continuum is to heal the surface wave field produced 
by compact sources, and therefore to hide them from view. It is found that the continuous spectrum 
is not a significant contributor to observable inter-ridge seismic power. 
Subject headings: Sun: helioseismology - Sun: oscillations 
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1. INTRODUCTION 

Acoustic oscillations in the solar convection zone are believed to be turbulently excited in a shallow sub-surface layer 



( Goldreich et al.|1994 Stein fc Nordlund|2001 Samadi et al.|2003| . In acoustic glor ies surrounding act ive regions, local 
helioseismic imaging reveals very compact bead-like emission in the 5-6 mHz band ( Donea et al.|2000" I , suggesting that 



at least a component of wave excitation is coming from effectively point-like sources with spatial extent far smaller 
than the typical wavelengths they generate. It is unclear whether these sources are associated with thin magnetic flux 
tubes, but such small scale structure does suggest the possibility. 

The thin m agnetic flux tubes that ubiquitously thread the solar surface are known to be potential acoustic sources 
(Cally 1986). Essentially transverse slow magneto-acoustic waves travelling down these tubes into the Sun's interior 
exhibit ever smaller wavelengths as the Alfven speed decreases with depth ( | Cally et al. 1994 Bogdan et al. 1996). 
Matching transverse velocity across the tube boundary therefore results in the generation of acoustic waves in the 
surrounding non-magnetic atmosphere with very large vertical wave numbers k z . The horizontal wavenumber kh, 
given by the acoustic dispersion relation fc 2 = (cj 2 — w 2 ) / 'c 2 — k\ (where uj c is the acoustic cutoff frequency) , is therefore 
necessarily imagi nary for given frequenc y ui and sound speed c such that fc 2 > (a; 2 — w 2 )/c 2 . This manifests as the 
'acoustic jacket' (Bogdan & Cally 19951, a horizontally evanescent sheath surrounding the flux tube. The acoustic 
jacket is distinct from the usual discrete spectrum of /- and p- modes familiar in helioseismology, and mathematically 
is associated with an additional continuous spectrum. 

The acoustic jacket is widely known but difficult to work with, though of course it appears naturally in grid- 
based simulations. Its mathematics is rather challenging, and not easily implemented in modally-based calculations. 
Indeed, it has become common practice to artificially truncate the atmosphere at some dept h so as to introduce a new 
unphysical discrete set of modes that can partially mimic the continuum of jacket modes (Barnes & Cally|[2000[ ), or 
even to ignore the acoustic jacket totally (with apologies) deeming it too difficult (e.g., Jain et al.||2(J(J9[ ). 

Compact sources are easily introduced into simulation codes, and acoustic jackets will appear naturally, though 
spatial resolution is an issue if they are thin. However, there is some value in understanding them from an analytic 
point of view to better alert us to their potential consequences. In that vein, the purpose of this paper is to revisit 
the basic jacket mode theory and show how vertical line sources may be constructed with it. This is then used to 
construct full 3D cylindrical or 2D cartesian wave fields originating from (specifically) point sources, which naturally 
exhibit a rich spectrum of high-fc z oscillations. We pay particular attention to the observable /- and p-mode surface 
power as functions of source depth and frequency. It is found that the continuum acts to 'heal' the surface signature 
of point sources by 'filling in' their image in the discrete-spectrum surface velocity field. 



EIGENFUNCTION EXPANSION OF LINEAR WAVES IN A TRUNCATED POLYTROPE 

2.1. Discrete and Continuous Spectra 
Details of the develo pment of the wave equations in a truncated adiabatic polytropic atmosphere may be found in 



Bogdan & Cally (1995 Section 2). We summarise them here. 
Consider an adiabatic plane polytrope of index m = 1/(7 



1), where 7 is the ratio of specific heats, such that the 
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equilibrium density and pressure vary with depth z according to 



P(z) = Pa 



P{z) = Pa 



m+1 



(i) 



where pa — pogzo/(m + 1), g is the gravitational acceleration, and zq is the truncation depth, i.e. zq z < oo. A 
stress-free boundary condition = is imposed at zq, where £ is the fluid displacement vector, and regularity is 
enforced at z — oo. 

This atmosphere supports both /- and p-modes, but not g-modes as it is adiabatic (the Brunt- Valsala, frequency 
N = 0). The vorticity therefore vanishes, and hence the displacement is potential, £ = ZqV$ (the factor of z 2 has 
been included to make $ dimensionless). Adopting the convenient non-dimensionalization s — z/z ^ 1, the wave 
equation takes the form 



m d 
s ds 



v 
s 



$(xi,s) = 



(2) 



where V = zoV and u = u> \J mzo / g is dimensionless frequency. The velocity potential $ may also be understood as 
a simple multiple of the Eulerian pressure perturbation, 



The boundary conditions are 



d_ 
ds 



v 
m 



$ = 0, as s ->• 1+, 



$ -> 0, as s 



(3) 



(4) 
(5) 



Equation ^ with boundary conditions (451 may be analysed using separation of variables. The vertical behaviour 
is found to satisfy the ordinary differential eigenvalue equation 



d^ 
ds 2 



m d 
s ds 



+ A $ 



where a monochromatic potential 



s, t) = $|| (s) exp[i(kx - ut)] 



(6) 



(7) 



is assumed, with x being horizontal distance and t time and A = — k 2 Zg. For simplicity we have assumed a two- 
dimensional [x-z] scenario; cylindrical 3D is easily addressed using Bessel functions in horizontal distance rather than 

a complex exponential. More general time dependence may be recovered by Fourier superposition. 

Equation Q with the stated boundary conditions may be solved in terms of Whittaker functions ( DLMF|2012 Ch. 
13). Specifically, we define 



(s;A) 



-l/2-H 



2 M M Ki _ M _ 1 (C)M K> „«) + 



(1/2 + /i) 2 - k s 



V( S ;A) 



2(l + /i)(l + 2^i) 



2 M K)M+ i 



(C)M K ,_ M ( S C) 



(8) 



(9) 



where fi = (m — l)/2, k — iv 2 /2\/A, and ^ = ~2iVX. In numerical examples, we shall use [i = 1/4 throughout, 
corresponding to m = 3/2, i.e. 7 = 5/3. It is assumed that k and £ are defined on ^ argA < 2-7T, so that both argK 

and argC lie in [— tt/2, 7r/2). (Computationally, it is convenient to express these as £ = 2/y — 1/A and k = ^ 2 /C where 
the square root applies the usual cut along the negative real axis.) These two solutions of Equation ([6]) satisfy the 
boundary conditions at s = 1 and 00 respectively. Their Wronskian 

(10) 



W A [0,V](s) = 2 l ^- 1 - 2 ^ (a/A - v%) W/,,m+i(0 



possesses simple zeros at \/Ao = iv 2 /{l + 2/i), corresponding to the /-mode w 2 = fc^, and an infinite sequence A n 
(n = 1,2,3,...) of p-modes such that Ao<Ai<A2 <•••<(). These become simple poles in the Green's function, 
a nd form the familiar dis crete spectrum. 

|Bogdan fc Cally ( 1995 1 go on to show that the /- and p-modes are not complete. One must also include a continuous 
spectrum, representing the jacket modes. Specifically, a suitably integrable arbitrary function f(s) may be expanded 
thus: 



71=0 



F n 1p n (s) 



T\ 4>(s; A) dX 



(11) 
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where ip n ( s ) = ^( s i A n ) an d 



r 1+a »Mr)f(r)dr, (n = 0, 1, 2, . . .), 



Here 



^A = -4a y r 1+2 ^(r; X)f(r) dr. 



-2i A A^ MKO '" M " l(Co) 



W«o,/i+l(Co) 
» r(l/2 - /i - K») M Kn _ M _i(Cn) 

^-7^ r(-2 M ) a/aA[VF K , M+1 (c)]U=A. 



(n = l,2,...) 



/A exp[i/ 2 7r/2v / A] 
4tt M 2 (A- A ) 



^+i(or 



(12) 
(13) 

(14) 
(15) 
(16) 



A,. 



with K n and £ n being the values at A 

To put our dimensionless units in a solar context, we adopt the numbers quoted by Bogdan & Cally (1995) for the 
m = I polytrope with po = 1.96 x 10~ 4 kgm -3 , po = 6.20 kPa, and g = 277.5 kms -2 : viz. zo = 285 km, from which 
a frequency of 1 mHz corresponds to v = 0.247 and 5 mHz to ^ = 1.23. Dimensionless frequency v = 1 corresponds to 
dimensional frequency 4.05 mHz. Frequencies are therefore seismically interesting, and sources placed at s = 2 

are only 285 km below the surface. For comparison, the detector resolution per pixel of the Helioseismic and Magnetic 
Imager (HMI) aboard the Solar Dynamics Observatory (SDO) is around 360 km. 

2.2. Special Cases: Orthogonality 

The expected orthogonality is easily confirmed. Specifically, with f(s) — ip n (s) we find that Fj = Sj n an d F\ = 0. 
Similarly, for f(s) = cf>(s; A) it follows that F n = for all n and F\ = S(X — A) (seelHindman & Jain||2012 Appendix 
A). 

2.3. Special Case: Delta Function 

Unlike the familiar case of an oscillating string of finite length, arbitrarily fine spatial scales are not attained by the 
discrete spectrum of waves in a polytrope. This is because, as n increases jthe eigenf unctions push deeper rather than 
becoming more oscillatory (per unit length). This is illustrated in Figure IT] It is therefore impossible to construct an 
arbitrary structure from the /- and p-modes alone, particularly a structure of small extent. 



1.0 F 




Figure 1. The /- (full curve) and first three p-modes (respectively long dashed, short dashed, and dot-dashed) as functions of dimensionless 
depth for the case fi = ^, v = 1. Higher order modes extend ever deeper rather than displaying finer structure. The modes have been 
normalized to (—1)™ at the surface. 



The most extreme example is the Dirac delta function. With f(s) = S(s — a), (a > 1), Equation (11) yields 



Coo ^oo \ 

J2A n M*)Mo-)+ J o Ax4>(s-X)4>(a-\)d\\ 



(17) 
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It is interesting to see how the continuous spectrum does this mathematically. Clearly the <5-function requires indefi- 
nitely large wave numbers to synthesise its sharp peak. This implicates the large A behaviour of A\cfr(s; X)cj>(a; A). It 
is easily shown that 



s;A) = ( 2fis-^ 2 -^ + O 



1 

7x 



cos 







(to: 









V\{s - 1) + o 



1 

7X 



(18) 



as A — > oo (using formulae from DLMF|2012 Section 13.19). In the same limit, A\ ~ I / (4tt /j, 2 y/~X) , so the dominant 
contribution to <j 1+2fl A\(f>(s; X)(j)((r; A) is 



"A 



1 /(7\l/ 2 +M ._ / _ 

(s, cr) = = I — COS V A(s — 1) cos V A(er — 1). 

7rVA Vs/ 



(19) 



This is readily integrated using a simple change of variables and the standard result J °° cos fcx cos /c£ dk = ^S(x — £) 
to yield 



<za(s, a) dX — S(s — a). 



(20) 



Hence 



/ [A x <j>{s; X)<f)(cr; A) - a-^a^s, a)] dX + Y] A n ^ n {s)^ n {a) - 0. 

J n=0 



(21) 



This makes explicit how the continuous spectrum forms a (5-function whilst also cancelli ng o ut the /- and p-modes. 

Although the integral on the left hand side is computationally problematic, Equation (21 1 is confirmed numerically 
to several significant figures. To avoid extreme round-off error and consequent catastrophic loss of precision it is best 
to integrate the two terms in square brackets separately. The integrals may be truncated, A < A max , at sufficiently 
large A max since the two terms are barely distinguishable there. A max = 5000 is generally adequate. The discr ete sum 



can also be very slowly converging in practice, but is easily accelerated using repeated Shanks transformation (Bender 
& Orszag 1978). In practice, using n = 0, . . . , 20 with ten Shanks transforms works as well as retaining 120 or more 
terms without transformation, and is much quicker to evaluate. 



2.4. Special Case: Gaussian 

We next consider a Gaussian g(s) = exp[— 20(s — 2) 2 ]. This is numerically simpler than the (5-function case, as the 
integral in Equation (11) converges in the normal sense. Figure [2] displays the reconstructed Gaussian, as well as the 
discrete and continuum contributions. Once again, as for the o-runction case, the sharp peak is contributed by the 
continuous spectrum. There is no hint of it in the discrete sum of /- and p-modes. 




Figure 2. The reconstructed Gaussian g(s) = exp[— 20(s — 2) ] (full curve, indistinguishable from the original function) together with the 
discrete (dotted, dominated by the /-mode) and continuum (dashed) contributions (see Equation 
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3. CONSTRUCTING A WAVE FIELD WITH A FORCING TERM 

3.1. Forced Wave Equation and Solutions 

The above mathematical development may be used to construct the full wave field <l>(x^, s) in multiple dimensions. 
Let x denote the horizontal coordinate in a 2D cartesian setting, and r the usual radial coordinate in an axially 
symmetric scenario. Both are assumed to have already been scaled by zq, and are therefore dimensionless. Consider 
a forcing term at x — or r = on the right hand side of the wave equation (assuming an e~ lVT time dependence as 

usual, where r = t>J g/mzo is dimensionless time), 



d_ 
~ds~ 



I f(s)S(x) (2D cartesian) 

*(xi,*) = /( s )a(xi)= . 

7r r 1 f{s)S(r) (cylindrical) 



(22) 



Only axisymmetric drivers are considered in the cylindrical case. The forced solutions may be written in terms of the 
discrete and continuous spectra: 



$f(x, s,t) 



1 -ivT 

2 e 



Jo v A 



dX 



or 



®t(r, s,t) = -e" 



— n 



(23) 



(24) 



where /Jq 1,1 = Jo + i Yq is the Hankel function of the first kind of order (representing an outward travelling wave) 
and K is the modified Bessel function of the second kind (representing the jacket). Both are logarithmically singular 
on r = 0. (The factors multiplying the sum and integral in Equation (24 1 may be derived by integrating Equation 
(22) over a disk of arbitrary radius and applying the div ergence theorem and several Bessel function identities.) The 
coefficients F n and J-\ are as given by Equations (12 13). Written in these forms, the continuum contribution is real 
and the discrete contribution is complex. For the axisymmetric case there is a logarithmic singularity at r = in both 
the continuum term and the real part of the discrete term. The imaginary part of the discrete term is not singular. 
An alternate construction, where $ or d^/ds is prescribed directly on the axis, is discussed briefly in Appendix [A] 
We now restrict attention to the axisymmetric point source, f(s) = 5(s — a), a > 1. The (5-function case is rather 
tricky numerically. It is again best to extract the large-A behaviour from the integrand and sum it exactly in the sense 
of generalized functions. To that end we rewrite the continuum integral (with T\ = o- 1+2 ^A\(j){o-; A)) in the form 



K {V\r)F X (j){s\\) dX 



1 



lim 



^/r 2 + (s + cr - 2) 2 Vr 2 + (s 
" K (VXr)A x 4>{o-, A) 0(s; A) dX - 




K Q (VXr)a\(s, a) dX 



(25) 



where the first term on the right hand side is the exact limit of the second integral in square brackets. It's form 
indicates that the dominant high-wavenumber contribution to the continuum integral is simply ra dia l inverse-distance 
falloff from the source at (0,cr) and its image at (0,2 — a). As before, a\ is given by Equation (19). Then both the 
first and second integrals need only be integrated numerically over < A < A max for sufficiently large but in practice 
moderate A max ((10/r) 2 is typically more than adequate), since their integrands nearly cancel at large A. This device 
significantly reduces numerical noise at small r, but is unnecessary at large r where the decaying Bessel- Kq accelerates 
convergence. 

Figure [3] and accompanying animations illustrate the generation of /- and pi modes by axisymmetric 5-function 
sources oftwo depths (and frequencies; the second panel corresponds to a higher frequency so as to make the pi mode 
shallower for graphical purposes). 



3.2. Surface Velocity 

From a seismic perspective, we may be most interested in the observed vertical component of velocity at the surface, 
v z = —ivd<fr/ds (where a dimensional scaling yj gzo / m has been removed). Since d<fr/ds = — (^ 2 /m)$ at the surface 
by boundary condition (HI), we may identify v z with $ there, apart from a normalization. 

The contour plots in Figure [3] and accompanying animations suggest that the /-mode is preferentially excited by 
shallow sources, and the p-modes by deeper excitation, presumably because the /-mode eigenfunction is dominant near 
the surface, with the p-modes residing considerably deeper (see Figure mi. This poses the question of which modes 
precisely are excited by a frequency v point source at depth a. At sufficiently large radius (beyond the jacket) the 
surface velocity power is dominated by the discrete spectrum and the Hankel functions become essentially complex 




Figure 3. Animations (contour plots of <J>f ) for the axisymmetric point source with v = 1, a = 2 (left) and u = 2, a = 4 (right), illustrating 
the generation of respectively the /-mode by a shallow source and the pi mode by a deeper source. Blue and red shading indicates different 
signs of <&. The colour maps and contour sequences are clipped near the source. Run with 'Loop' set. Files: Deltad_l_2.mp4 and 
Deltad_2_4.mp4. 

exponentials (with an r -1 / 2 amplitude factor), so we have the usual Fourier orthogonality in the horizontal plane, 
allowing us to write velocity power atQ 



1 00 1 
P(r)oc-£— L = \F n ^ n {l)Y 



(26) 



This may be summed, and the fractions in each of the modes n, with n = corresponding to the /-mode and n > 
to p n , may be calculated for various frequencies and source depths. The results are summarised in Figure [4] 
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Figure 4. Fractional surface power in the / (full curve), pi (long dash), p2 (short dash), p3 (dot-dashed), and p4 (dotted) modes excited 
by a harmonic axisymmetric <5-function source of frequency v placed at depth a. Left: u = lj right: v = 2. In each graph, the cr is 
arbitrarily chosen to range down to the p4 Lamb depth, = —v 2 /\^. We see that fractional power is almost exclusively dependent on n 
and cr/sr, alone. 



1 The I-Hq 1 ^ (k n r)| 2 term decays with radius as 27r 1 fc n 1 r 1 for sufficiently large r . 
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The dependence on frequency may be understood by recourse to the complete polytrope, which extends all the way 
up to z = and is a pure power law with no length scale. Its discrete eigenvalues are exactly w 2 = gk{\ + 2n/m). The 
depth dependence of its eigenf unctions therefore scales exactly with k~ l and hence with oj . Consequently, p-mode 
excitation as a function of point source depth can depend only on depth through kz or equivalently urz. For the 
truncated polytrope this is only approximately true. Nevertheless, Figure [4] indicates that it is a good approximation. 

Figure [5] depicts surface v z as a function of radius for four source depths of frequency v = 1. It is clear that the 
acoustic jacket at the surface extends ever farther as the source moves deeper. In each case there is of course a very 
compact logarithmic singularity in both the discrete and continuum velocities, but these cancel exactly in the full 
solution, since the current case is equivalent to specifying 27r$(0, s, t)/ lnr — > f{s)e~' lVT as r — > (see Appendix |A[), 
which vanishes at the surface for f(s) = 5(s — a). This does not mean that $ vanishes at r = 0, s = 1, only that the 
singular part does. That leaves only nonsingular components. We refer to this as 'healing' the logarithmic singularity, 
though the a = 1.2 case in Figure [5] illustrates that the healing can still leave a sharp velocity peak at r — if the 
source is too shallow, as one might expect. The removal of the logarithmic singularity still leaves a broader continuum 
jacket as is apparent from the grey curve in the figure. Apart from the very shallow source case, it is interesting that 
even this 'outer' jacket seems to largely counteract a similar feat ure in the real part of the discrete solution, further 
hiding the source from view. In the 2D cartesian case (Equation (23)), where there is no singularity, only this regular 
part of the jacket exists. 
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Figure 5. Surface vertical velocity v z against radius r for 5-function sources at four different depths: a = 1.2 (top left); <r = 2 (top right); 
a = 6 (bottom left); and a = 12 (bottom right). Note that the a = 1.2 case is depicted over a reduced radius to better display its behaviour 
near the origin. In each case the thick black full and dotted curves represent real and imaginary parts of the full solution. The grey curve 
represents the continuum velocity (which is purely real). The dashed curve is the real part of the discrete component (the imaginary part 
is again the thick dotted curve). 



The dependence of surface signature on frequency is illustrated in Figure [6] for a shallow point source at a — 1.2, 
i.e. only some 57 km below the surface. At high frequency v — 2 (8.1 mHz) the discrete spectrum dominates since 
the jacket is extremely narrow (grey curve), and merely cancels the singularity in the discrete solution (dashed). On 
the other hand, small-r surface velocities are progressively dominated by the continuous spectrum as v diminishes. 
Since the dominant behaviour a a of that contribution is independent of v (see Equation (19)), we expect only weak 
dependence on frequency of v z at small r as v — > 0. This is illustrated by the v — 0.5 (2 mHzJand 0.25 (1 mHz) cases in 
Figure [6] being very similar. In both instances, there is a prominent and compact surface velocity feature immediately 
above the source. It is produced by the continuum (grey curve), suggesting that shallow low frequency jackets are 
potentially visible. The near-field is made more apparent at low frequencies because of the very low amplitude of the 
far-field. As the source is moved deeper, the ability of the continuum to heal and hide comes to the fore, and the 
surface 'bump' becomes much less prominent. 

Despite the acoustic jacket being composed of a continuum of horizontal wavenumbers, its effect at the surface is to 
heal the logarithmic singularity and make further relatively minor changes to v z within a small radius about the source 
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Figure 6. As for the a = 1.2 case of Figure [5[ but with frequencies v = 2 (top left), v = 1 (top right), v = 0.5 (bottom left) and v = 0.25 
(bottom right). In each case the thick black full and dotted curves represent real and imaginary parts of the full solution respectively. The 
grey curve represents the continuum velocity. The dashed curve is the real part of the discrete component. 



axis. Nevertheless, there may still be some measure of continuum power observable that shows up in the seismic power 
maps. To look for this we run a 'slit' across the disk through the axis of the source, and compute the vertical velocity 
power spectrum for two source depths (Figure [7]). Using HMI-like per pixel spatial resolution Ar = 1.27, there is no 
sign of any significant inter-ridge power in either case. The sources have very effectively hidden themselves. 
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Figure 7. Surface velocity power spectra of \/r 2 + 1 v z along a slit through the source axes for v = 1 sources at a = 2 (left) and a = 6 
(right). The factor \/r 2 + 1 removes the l/y/r geometric amplitude attenuation without killing off the signal at r = 0. We use a spatial 
grid of 1024 points of spacing Ar = 1.27 (comparable to HMI resolution) and apply Hann windowing to suppress spurious side lobes. The 
vertical grey lines represent the eigen-wavenumbers of the /- and first five p- modes (right to left). As we saw earlier, the shallow source 
overwhelmingly excites the /-mode (peak power around 2.4 at this resolution), whereas the deeper source excites substantial / and pi 
power, and even noticeable p2, ■ ■ - P4- No significant background power is apparent. The /-mode peak is cut off for graphical purposes in 
the left panel, but extends up to about 2.4 at this resolution. The normalization of the two sources is the same. 



4. CONCLUSIONS 

The acoustic continuum is an intrinsic component of the Sun's interior wave field, partic ularly for small scale sou rces 
not resolvable by summation over the discrete modes. We have summarized the theory of Bogdan & Cally ( 1995 1 and 
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applied it to constructing compact sources, in particular point sources. This has allowed us to visualise the role of the 
near-field in generating /- and p-modes, and to quantify the relative proportions of surface velocity power going into 
each of the discrete modes from monochromatic point sources of varying frequency and depth. As expected, shallow 
sources preferentially excite the /-mode, with subsequent p-modes taking precedence as the source is moved deeper or 
the frequency increased. 

The significance of the continuum 'acoustic jacket' for observationally accessible surface velocities is explored in 
Section |3.2| Within a small radius determined by source depth, the continuum plays a significant role. It heals the 
coordinate singularity at r = 0, but does so with (of course) a continuum of wavenumbers. (In the real Sun, compact 
acoustic sources will have some non-zero spatial extent, possibly hundreds of km or more, making the logarithmic 
singularity irrelevant, but the 'outer' jacket will persist.) The resulting surface v z displays a simple sinusoidal behaviour 
that is basically unremarkable, except if the source is very shallow indeed (see the a = 1.2 case in Figure [5] and the 
low frequency cases of Figure pi. The effects of the continuous spectrum for buried point sources do not contribute 
substantially to background power in the surface doppler power spectrum for the cases considered. 

The fact that the wave field close to the source beneath the surface is dominat ed by the jacket modes must surely 
place limits on how well reconstructive techniques such as Acoustic Holography (Lindsey & Braun 2000) can image 
them. We have not attempted this, but the optical assumptions underlying the method clearly break down within the 
jacket, irrespective of any further finite wavelength effects. 



APPENDIX 

WAVE FIELD MATCHING TO A PRESCRIBED SAUSAGE OR KINK OSCILLATION 

An alternate scenario relevant to a thin vertical magnetic flux tube at r = is to directly prescribe the pressure 
perturbation $ (the axisymmetric sausage mode) or the radial displacement <9$/<9r (the kink wave, with an e ±%e 
azimuthal dependence) on the axis. For the kink, with 2nr 2 d<S>/dr -> f(s) e^ 9 '"^ as r -> 0, 



$ k M,s,r) 



V^tff 3 (^Kr) F„V„(s) 



n=0 



If 00 f \ 



(Al) 



where we have used H^\{z) = —H^ l \z) and K-i{z) = K\{z). For the sausage mode with 27r$(0, s, r)/ lnr 
f(s) e-' lVT as r -> 0, 



$ s (r,s,r) 



Y, #o 1} (^Kr) F n ^ n {s) + — K (VaY) T x <f>{s; X) dX 

m — n J 



(A2) 



$ f (r, s,t) ■ 



We see that the axisymmetric forced case and the sausage case are equivalent. We will not pursue the kink case here, 
except to note that the added factor of vA in the integrand emphasises higher wavenumbers, meaning that the kink 
jacket is more tightly bound to the source. 



The author thanks Sergey Zharkov for asking the question. 
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